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. The Oslo rice pile model is a sandpile-like paradigmatic model of "Self-Organized Criticality" 

' (SOC). In this paper it is shown that the Oslo model is in fact exactly a discrete realization of 

the much studied quenched Edwards- Wilkinson equation (qEW) [Nattermann et al., J. Phys. II 
France 2, 1483 (1992)]. This is possible by choosing the correct dynamical variable and identifying 
' its equation of motion. It establishes for the first time an exact link between SOC models and the 

I field of interface growth with quenched disorder. This connection is obviously very encouraging as 

, it suggests that established theoretical techniques can be brought to bear with full strength on some 

' of the hitherto elusive problems of SOC. 

(N : 

PACS numbers: 64.60.Ht, 05.65.-|-b, 68.35.Fx, 02.50.-r 

O ' The Oslo rice pile model (Oslo model hereafter) was originally intended to model the relaxation processes in real 

, rice piles l]|. Meanwhile, it has been subject to many investigations and publications in its own right. The model as 
defined below supposedly develops into a scale free state without the explicit tuning of external parameters, and is 

, therefore regarded as an example of Self-Organized Criticality (SOC) In fact, contrary to many other "standard" 

^ ; models of SOC 0, i, S , it shows a reliable and consistent (simple) scaling behavior and is robust against certain 

^ , changes in the details of the dynamics |2, i^iSj- The most prominent observable in the model, the avalanche size s, is 

' governed by a probability distribution V{s) which obeys simple scaling, 

S : Vis) = s-^g{s/so) and so = , (1) 

(— ( where L denotes the system size and r and D are critical exponents, consistently reported to be r = 1.55(10) and 
Q D ^ 2.25(10) 0, S H O m E3- These two exponents are related hy D{2 - t) = 1 ^E^, which can be proven 
O , easily given that the first moment of 'P(s), (s), scales like L. 

In the following the model is defined, the relevant dynamical variable extracted and its equation of motion derived, 
fSJ , which turns out to be a discretized quenched Edwards- Wilkinson (qEW) equation. By analyzing the essential charac- 
^ ■ teristics of the model on the lattice, such as uniqueness of the solution and symmetries, it is then possible to construct 
[ the continuum theory, which can subsequently be examined using standard methods. 

■ The model 0| is defined on a one dimensional grid of size L, where each site i — 1 ■ ■ ■ L has slope Zi and critical 
Qv^ , slope zf £ {1,2}. Starting from an initial configuration with Zi = and random everywhere, the model evolves 

■ according to the following update rules: 1) (Driving) Increase zi by one. 2) (Toppling) If there is an i with Zi > zf 
OA , decrease Zi by 2 and increase its nearest neighbors by one, Zi±i — > Zi±i + 1, provided that 1 < i ± 1 < L. A new 

' zf is chosen at random, 1 with probability p and 2 with probability q = 1 — p. 3) Repeat the second step until 
Zi < zf everywhere. Then proceed with the first step. The order of updates is irrelevant in this model and the original 
definition does not fix it explicitly. Therefore the microscopic (fast) timescale is a priori undefined. 
^ , The avalanche size s is defined as the number of charges, i.e. apart from boundary effects, twice the number of 
• times the second rule is applied between two consecutive application of the first rule. For convenience the model is 
^ [ dissipative on both boundaries, where one of the two "units" lost by the boundary site during toppling leaves the 
O ■ system. 

A few years ago Paczuski and Boettcher translated the Oslo model into the language of interfaces in random media 
However, the evolution of the dynamical variable H{x,t), which is the total number of topplings of site x, was 
given by dtH = 0{d'^H — ?]{x, H)), where dt is defined in discrete time, i.e. dtH = H{x, t -|- 1) — H{x, t) and is the 
5^ • lattice Laplacian, so that x is actually an index. The last term rj^x, H) represents a quenched noise. The Heaviside 
\ 0-function makes this equation of motion highly nonlinear and analytically almost intractable |ia|. Paczuski and 
Boettcher have already conjectured that the Oslo model is in the same universality class as qEW [l^. More recently, 
Alava has suggested that certain other sandpile models are described by qEW 0. It is, however, important to realize 
that no rigorous and exact link has so far been established between SOC models and the qEW equation. 

The crucial step to make this correspondence exact is to identify the proper dynamical variable. It is found in the 
form of the number of times a site has been charged (i.e. received a unit from a neighbor during a toppling or by 
external drive, see below) h{x,t), where x and t are discrete for the time being. There is a simple functional relation 
between h{x,t) and H{x,t), which can be obtained as follows: Each site can be in one of three stable configurations, 
Zi S 0, 1, 2. When a site receives a unit from a neighbor, it changes state as shown in Fig.^ Charging a site in state 



FIG. 1: Each site can be in one of three states and changes stepwise between them, whenever it receives a charge. The labels 
indicate the probability of the move and whether it entails a toppling. 



necessarily leads to state 1 without toppling and the specific value of z'i is completely irrelevant at this stage. Similar 
for state 2: If a site receives a charge in this state, its must be 2 and it must topple. The only point where the 
value of z'i actually matters, is in state 1, therefore it can be effectively chosen at random when necessary, so that 
the site topples with probability p (according to the probability of having = 1) or increases to 2 with probability q 
(see Fig.^l. It is immediately clear that any even number of charges, say m = 2n, starting from Zi = \ leads to state 
1 again with n topplings. An odd number of charges, say m = 2n + 1, leads either to n topplings and state 2 or n + 1 
topplings and state 0. This is illustrated in Fig. The m charges lead to m steps along the arrows. Whenever one 
moves left, the site topples. 

In order to write a functional relation between h{x, t) and H(x^ t), the randomness in the decision of moving to the 
left or to the right from state 1 must be quenched in h{x,t), i.e. it is not allowed to change unless h{x,t) changes. 
This can be summarized as 

H{x,t+1) = ^{h{x,t)+r]{x,h{x,t))) , (2) 

where 77 is whenever h{x,t) is even, corresponding to state 1. If h{x,t) is odd, 77 is either 1 (with probability p, 
state = 0) or — 1 (z, = 2). Every sequence of ri{x,h) values maps uniquely to a sequence of z^ and vice versa. 
The equation above can easily be transformed to comply to any initial configuration, especially to Zi{t = 0) = 0. 
Essentially, it is (O, which makes the exact identification of the Oslo model and qEW possible. 

The final equation is derived by noting that obviously /i(x, t) = H{x — 1, i) + H{x + 1, t) with appropriately chosen 
boundary conditions (BC's) (see below), so that using the short hand notation — h[x± 1, t) and i]^ = r;(a;± 1, h^) 
the equation of motion is 

h{x,t+l) -h{x,t) = ]^{h~-2h{x,t) + h+ + i]++'q~) , (3) 

which is the exact representation of the Oslo model as defined above, captured in a single equation. Its differential 
form is accordingly 

dth{x, t) = ^dlhix, t)+(^l + ijix, hix, t)) . (4) 

The right hand BC is h{x — L + l,t) = (and h{x — L,t) = in the continuum), while the left hand EC provides 
the driving via h{x = 0,t) = 2E{t), E{t) being the total number of initial seeds (step 1 above) at time t. These seeds 
arrive at site x = 1 via the Laplacian. In the continuum, the simplest drive is E(t) = vt with v a driving velocity and 
t the microscopic time. Together with the BC's, Eq. Q or the generalized form 

dth{x, t) = Kdlh{x, t) + .g ^1 + vi^i H^, t)) , (5) 

where the correlator of 77 is now normalized, i.e. J dx J dh{riri) = 1, describes the movement of an elastic band over 
a rough surface '16*1 pulled by a transverse force acting at one end point only. Below it is shown that the A-term 
disappears in the continuum, establishing the first rigorous identification of the Oslo model and the qEW equation. 
The same equation with different properties of the noise term and/or different BC's applies to other models, such 
as the BTW model Fixed Energy Sandpiles (for example ^3) or the tilted sandpile Having identified the 
relevant dynamical variable h, the effect of modifications of the dynamical rules of the Oslo model, such as |[2lSi3]j 
can be understood. 

The equation above exemplifies a general "trick" |29| to get rid of ^-functions in equations of motion — they often 
appear in descriptions of sandpile-like systems (for example One simply replaces 9{h — he) hy h + rj{h) with 
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an appropriately chosen sawtooth-like 77. This does not necessarily simplify the problem, unless there is already a 
quenched noise present in the system. In this case the 9 turns into a correlation in 77. This is highly remarkable from 
the point of view of SOC, because the presence of "thresholds" is usually expected to be a crucial ingredient of SOC 
HiQiH^' Moreover, the correlations in 77, which are of fundamental significance in interface models 0,1^ and 
have been neglected in former mappings, now arise naturally from the dynamical description of the model. 

In order to construct the proper continuum theory, it is worthwhile to consider the formal solution of Eq. I^l. It 
will turn out later that E{t) = vt is sufficiently general, so that it makes sense to define v{x) = v^^^ and 

h{x,t) = 2v{x)t + P3{x) + z{x,t) (6) 

in order to homogenize the BC's. ^3(2;) is a third order polynomial only present to cancel the first term in the 
differential equation, i.e. k&^P^ — 2v{x), with roots at x = and x = L. Therefore dtz = nd^z + grix(x,h{x,t)) 
with homogenous BC's. The term ri\{x,h{x,t))) = (1 + \^^)ri{x,h{x,t)) is actually a functional of h. The initial 
condition of z{x, t) is not z(a;, t = 0) = as for h, because of the data shift above. But due to the homogenous BC's 
any initial condition decays, so that the initial sources, accounting for z{x,t = 0) = —P^{x)^ can be ignored. Then 
the formal solution is z{x,t) = ^n(i) sin(A;„a:) with 

Zn{t) = T" / dt' ( dx'r]x{x',2v{x)t + z{x\t)) 
^ Jo Jo 

X sm{knx) exp(— A;^k(< — t')) (7) 

and kn = ^. 

According to Eq. the tilt of h[x,t) in x increases in time. Assuming stationarity of the relevant statistical 
properties (especially avalanches as defined below), this requires the solution to be invariant under tilt, which is also 
known as Galilean invariance |2^ : h' = h + ax must produce the same statistics as h, which entails ri{x,a + ax) 
to be equally likely as ri{x,a), so that {'r]{x,a + ax)ri[x' ,a' + ax') — {r](x,a)r](x' ,a')). But assuming the standard 
form 14] {rj{x,a)rj{x' ,a')) = A||(a; — a;')Ax(a — a'), the correlator obeys for any x — x' where A||(x — x') is finite, 
A_L(a — a') = A_L(a — a' + a{x — x')). This holds for any a, so if A|| {x — x') was finite for any x — x' 7^ 0, A^ would 
be bound to be a constant. This is impossible, because A_l must be non-vanishing somewhere and normalizable, so 
that A||(a; — x') must vanish for any finite x — cc', i.e. it must be a (5-function. 

Next it can be shown that the Oslo model obeys Middleton's no-passing 'is']. For A 7^ this will lead to a 
constraint on the noise which is incompatible with the 5 correlation of A|| in the continuum, so that A must vanish 
in the contimmm. Defining a partial ordering >z_ for two configurations hi{ti,x) and h2{t2,x) of the interfaces as 
>z h2{t2,x) <^ ^ x&[o,L]hi{ti,x) > h2{t2,x), ouc has to show that this order is preserved under the dynamics 
|24j . With the "external field" being the BC's Ei{t) and E2{t), one shows that if hi{to,x) >: h2{to,x) for a given to 
(which entails Ei{to) > £^2(^0)) the interfaces can never "overtake" each other a,t t > to. By assuming the opposite, 
one only needs to prove that where the two interfaces "touch" for the first time, xo, the velocity of hi is higher or 
equal to the velocity of /i2- For the model on the lattice this is equivalent to 

h+ + 4 + + r/f >h++ r/2+ + + (8) 

using the same notation as in 0. In the original discrete model, condition (8) follows immediately from ri{x, h) + h 
being a monotonically increasing function in h for any x. For the continuum equation 10) the corresponding calculation 
gives 

\gdh-q{x,h)>~K (9) 

assuming that = dl-q + dxhdxdhrj + dxhdhdxr] -I- d^hdhrj + {dxh)'^df^r] and that the interface is smooth in xo such 
that dxhi{xo,t) = dxh2{xo,t) and d'^hi{xo,t) > dlh2{xo,t). For a noise with divergent width, A||(x) — S{x), Eq. @ 
cannot hold for any A 7^ 0, i.e. a non-vanishing A destroys no-passing. However, no-passing must be regarded as a 
crucial feature, as it ensures the asymptotic uniqueness of the configuration and is reminiscent of the irrelevance of 
the order of updates in the original model, so that A = is a necessary condition for the equivalence of the continuum 
and discrete model. 

This is physically justified: Assuming a smooth ij, in the continuum approximation of Eq. (jSJ A becomes proportional 
to the square of the lattice spacing and therefore vanishes in the continuum limit. 

Keeping the A term nevertheless, a nai've scaling analysis shows that it is irrelavant. Moreover, its Fourier transform 
in Eq. Q produces only a term —gXk^, because of the total derivative in -qx- This can be absorbed into the bare 
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propagator of a perturbative expansion in the style of [13 . l2ll | in the form '^^|^f.2 , leading possibly to an ultraviolet 
divergence. Apart from that, the terms obtained for an renormalization group treatment are structurally the same 
as in 21] as calculations show (details to be published later). The only differences are due to the peculiar way of 
driving the interface (i.e. the term 2v{x), which is a mean velocity in ((HJ, but also drives the model by moving the 
quenched noise in lO) and the non-conservative nature of the interface (which makes sense only for a finite system) 
leading to the homogenous BC's and therefore to the sin(fc„a;) rather than exp(2ifc„x) terms. In turn, the standard 
qEW problem [l4| corresponds to an Oslo model with periodic BC's and continuous, uniform drive. 
Expanding 77 in powers of Zn, the first two terms of Zn{uj) (the Fourier transform of |(7J| in t) are: 

2.9(1 - Afc^) f f\ ,.f , u; \ sin(/c„a;') 



L{nkl +iuj)\Jo 'V 2u(x') J 2v{x') 

[ dx [ dqy^ fl{x', q)^-^^—^t^—-^z„i{LU - 2v{x')q) sin(fc„x')) 
Jo J-00 „,_i v27r ^ 



m — 1 

where f]{x, q) is the Fourier transform of 77(2;, h) in h. 

The definition of the avalanche size s in the continuum is the area between the interface configurations at two times 
ti and t2, s — dx{h{x,t2) — h{x,ti)), so that (s) = vAtL with A< = t2 — ti, because {z[x,t)) is expected to be 
asymptotically independent of i, as a non-vanishing linit^oo d{z{x,t)) with homogenous BC's would require support 
for a divergent curvature of the interface. Choosing Ah = Atv constant for different system sizes L then preserves 
the property (s) cx L. 

Due to the asymptotic uniqueness of the solution the system can either be driven in jumps of Ah separated by 
sufficiently long times, or driven very slowly taking "snapshots" of the configuration in order to calculate s. 

The model possesses two characteristic timescales: One is the diffusive timescale to = L'^ /k, the other one is the 
non-trivial scale due to noise and drive, tg = /(v'^L). One has to maintain a sufficiently large At to prevent distinct 
avalanches from merging, otherwise the central limit theorem would turn V{s) into a Gaussian. The SOC limit is 
usually identified with w ^ 0, which makes sense only in the presence of an intrinsic scale for v. The only combination 
of parameters (k, g and L, but A — 0) which provides a "natural velocity" is Vg EE {g'^nf/^lL. The SOC condition 
?; ^ is therefore already met by w <C c>c , which is however, not sufficient. According to Ref. 11] At ^ 
with z K, 1.42, so that Aft, — const, entails v <C i^^, which therefore seems to be the correct condition for SOC, even 
though the microscopic timestep in 'llj is defined as a parallel update, which is not exactly JH)). 

Preliminary numerical studies indeed suggest that ((SJ with A = is a valid continuous description of the Oslo model: 
Fig. El compares a scaling collapse for different system sizes of the continuous model (with A = 0) and the original, 
discrete one. The best collapse is obtained by r = 1.55 for both models. The scaling law D = 1 + x [ill remains 
applicable as long as the two confi gur ations at ti and t2 are correlated. It is in perfect agreement with numerical 
results mill for the qEW model [H^. 

In conclusion, the Oslo model has been reduced to a quenched Edwards- Wilkinson equation. In the continuum 
limit the qEW becomes the exact equation of motion for the Oslo model. This not only makes it possible to approach 
the exponents of an SOC-model analytically, but also gives insight into the nature of avalanche like behavior and 
the relation between SOC and other theories of critical phenomena. It provides the perfect test bed for analytical 
methods proposed for SOC. 

The established relationship is presently being pursued in order to develop a direct approach to the critical exponent 
T, clear up the role of the noise and clarify the relation between noise and drive. The framework used here is also 
promising for other models, such as the BTW model 0], various other sandpile models (l&.lQj and the Zhang model 
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